SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

THEORY AND 
COMPUTATION 

METAL-ORGANIC FRAMEWORKS 



Received 
26 November 201 3 

Accepted 
28 January 201 4 

Published 
14 February 2014 



Correspondence and 
requests for materials 
should be addressed to 
Q.S. (sunqiang@pku. 

edu.cn) 



The superior catalytic CO oxidation 
capacity of a Cr-phthalocyanine porous 
sheet 



Yawei Li 1 &Qiang Sun 1 



1 Department of Materials Science and Engineering, College of Engineering, Peking University, Beijing 100871 , China, 2 Centerfor 
Applied Physics and Technology, Peking University, Beijing 100871, China. 

Two-dimensional organometallic sheets containing regularly and separately distributed transition atoms 
(TMs) have received tremendous attentions due to their flexibility in synthesis, well-defined geometry and 
the promising applications in hydrogen storage, electronic circuits, quantum Hall effect, and spintronics. 
Here for the first time we present a study on the superior catalytic CO oxidation capacity of a 
Cr-phthalocyanine porous sheet proceeding first via Langmuir-Hinshelwood (LH) mechanism and then via 
Eley-Rideal (ER) mechanism. Compared to the noble metal based catalysts or graphene supported catalysts, 
our studied system has following unique features: without poisoning effect and clustering problem, having 
comparable reaction energy barrier for low- temperature oxidation, and low cost for large-scale catalytic CO 
oxidation in industry. 

The catalytic oxidation of carbon monoxide (CO) has attracted extensive attentions for several decades due to 
its crucial role in exhaust treatment of the post-combustion process for automobiles 1 as well as in alleviating 
the poisoning effect on methanol fuel cell electrocatalysts 2 , it is also a prototype catalytic reaction in surface 
chemistry because of its stoichiometric simplicity as well as industrial significance, which plays an essential role in 
evaluating catalyst activity, selectivity and durability. Following the pioneering work on Pt wires of Langmuir 3 , 
many noble metals such as Pt 4 ' 5 , Rh 6 ' 7 , Pd 7 and Au nanoparticles 8 " 10 are studied for catalyzing CO oxidation. 
However, these catalysts are generally costly, and platinum based catalysts usually encounter CO poisoning 
problems, impeding durable operation. To overcome this problem, novel low dimensional nanocatalysts includ- 
ing metal embedded graphene 11 " 16 as well as metal embedded carbon nanotubes (CNTs) 16 have been suggested for 
enhancing the catalytic activity of CO oxidation with considerably low energy barriers 12 " 16 . However, it is 
technically difficult to controllably generate such defects in these materials. In fact, previous calculations revealed 
that the formation energy of a single vacancy (SV) or a double vacancy (DV) in graphene can be as high as about 
7.5 ~ 8.0 eV 17 . So it is highly desirable to develop novel porous organometallic materials intrinsically containing 
exposed metal cations for large-scale catalytic CO oxidation in industry. 

The recent success in embedding Fe atoms into a 2D polymeric phthalocyanine (Pc) sheet 18 inspired us to 
employ metal-phthalocyanine-based sheets as the candidates of low-temperature CO oxidation catalysts. Indeed, 
various experimental 19 " 24 and theoretical 25 " 27 efforts have been devoted to unveil the observed catalytic activities of 
molecular transition metal phthalocyanine compounds towards oxygen reduction reaction and C0 2 reduction 
reaction. The 2D Pc sheets, as illustrated in Figure 1, are superior to defected graphene/CNTs based systems due 
to the following reasons: First, the synthesis procedure of 2D FePc framework is flexible for other TM elements, 
indicating that other 2D TMPc sheets are able to be synthesized in the future. Second, the structural motif of 2D 
FePc framework highly resembles biological metalloporphyrin systems such as hemoglobin, where the kernel of 
hemoglobin possesses strong binding affinity to CO and 0 2 28 . Therefore, it seems plausible to assume that certain 
kinds of 2D TMPc sheets can also reversibly bind CO and 0 2 ligands, which acts as an indispensible prerequisite 
in the whole catalytic CO oxidation process. Third, the electronic states of TM atoms can be effectively tuned by 
the Pc square planar ligand field and by varying the metal atom at the Pc center 29 , enabling us to systematically 
investigate the catalytic behavior of various 2D TMPc sheets. Although the basic magnetic properties, intrinsic 
magnetic couplings and gas storage properties of 2D TMPc frameworks have been thoroughly investigated 30-35 , 
the possibility of such systems for catalyzing CO oxidation and their corresponding reaction mechanisms have 
not been well addressed yet. 
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Figure 1 | Schematic illustration of 2D TMPc sheet and its catalytic behavior toward CO oxidation. The highlighted areas denote the most catalytically 
active sites. 



In this work, we performed first-principles calculations to system- 
atically investigate the catalytic feature for CO oxidation of the 2D 
TMPc sheets. We chose TM atoms in the 3d series ranging from Cr to 
Cu. Among them, the 2D CrPc sheet is identified as a potential 
candidate for low-temperature CO oxidation, exhibiting the highest 
catalytic activity without the CO poisoning effect. Previous density 
functional calculations reveal that TiO 2 (110) or MgO(lOO) sup- 
ported model Au nanoparticle catalysts can accelerate CO low- tem- 
perature oxidation with energy barriers of 0.35 ~ 0.75 e V 9,36 ' 37 . For 
free-standing Au nanoparticles the energy barriers are even larger 
(0.46 ~ 1.03 eV) 38 . Using the same level of theory, we found that 2D 
CrPc can catalyze CO oxidation with an energy barrier of 0.55 eV, 
comparable with those of supported Au nanoparticles. The different 
extents of the hybridizations between CO/0 2 frontier orbitals and the 
TM d states are predominantly responsible for the observed trends in 
the adsorption energies of both CO and 0 2 on different 2D TMPc 
sheets. The reaction of CO and 0 2 on 2D CrPc sheets proceeds via 
two steps, namely the LH mechanism followed by the ER mech- 
anism. We further provided detailed analyses on the electronic struc- 
ture of the 2D CrPc system in the LH reaction processes to elucidate 
the superb catalytic activity exhibited in the 2D CrPc sheet. 

Results 

Searching for the candidate. Before examining the CO oxidation 
reaction paths on 2D TMPc sheets, we first studied CO and 0 2 
adsorption on these frameworks. The adsorption energy of an 
adsorbate is defined as £ ad = -Etmpc + moi — ^tmpc — E moh where 
^tmpc + mob ^tmpc an d £ mo i are me total energies of the 2D TMPc 
with the adsorbed molecule, the isolated 2D TMPc and the isolated 
gas molecule, respectively. By this definition, negative (positive) 
values of £ ad represent exothermic (endothermic) adsorption 
processes. The variation of £ ad (CO) and £ ad (0 2 ) in different 2D 
TMPc framework is illustrated in Figure 2. Both £ ad (CO) and 
£ a d(02) undergo a giant reduction in their absolute values when 
changing the central TM atom from Fe to Co. However, CO 
adsorption is appreciably enhanced while the binding to 0 2 
adsorption is only weakly attenuated when the TM atom is 
changed from Cr to Fe, resulting favored bindings to CO in MnPc 
and FePc 2D frameworks. 

From the above analyses we can see that among all of the studied 
systems, only the 2D CrPc, MnPc and FePc systems are likely to 
become the candidates of low-temperature CO oxidation catalysts, 
arising from the strong interaction between such 2D porous sheets 
and CO as well as 0 2 . Next we focus on another important factor 
which affects the durability of the CO oxidation catalysts, namely the 
CO poisoning effect. Because £ ad (CO) is larger than £ ad (0 2 ), the 



MnPc and FePc sheets will be dominantly covered by CO if the 
1 : 1 CO/0 2 mixture is injected with a constant gas flow rate. In this 
condition, the MnPc and FePc sheets adsorbed with CO will not 
suffer from the poisoning effect only when it is exothermic for those 
CO -covered 2D frameworks to bind 0 2 molecules and form coad- 
sorption states. To this end, we define £ coa d =£tmp c +co+o 2 — 
£tmpc+co-£o 2 > where £tmp c +co+o 2 > £tmp c + co and E 0l denote 
the total energies of the CO/0 2 coadsorbed on 2D TMPc, CO 
adsorbed 2D TMPc and free-standing 0 2 molecule, respectively. 
The negative values of £ coad would suggest the feasibility of removing 
CO from the 2D TMPc sheets. Our calculation results indicate that 
the FePc sheet is unable to form such a stable coadsorption state. 
Furthermore, although the coadsorption state can exist in the MnPc 
sheet, the £ coa d value of the MnPc sheet is 0.92 eV. Therefore, both of 
the two sheets are susceptible to CO poisoning in catalyzing the CO 
oxidation reaction, indicating that they cannot be employed as CO 
oxidation catalysts. In contrast to these two sheets, the catalytic active 
sites in a the CrPc sheet are much more stable towards irreversible 
CO binding with a negative £ coad of —0.16 eV, suggesting that only 
the CrPc sheet can be applied to catalytic oxidation of CO with 
negligible CO poisoning effect. 

The distinct adsorption behaviors of CO and 0 2 can be rationa- 
lized through the analysis of the spin -polarized orbital of 2D TMPc 
sheets at the Y point. The selected orbital energy levels of CrPc, FePc 
with their corresponding electron density distributions are illu- 
strated in Figure 3a and 3b. Generally, a chemical reaction cannot 
change the overall spin of the system due to the Wigner- Witmer rules 




Cr Mn Fe Co Ni Cu 

Figure 2 | Adsorption energies (E ad ) of CO and 0 2 adsorbed on 2D 
TMPc. 
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(a) 




(b) 



Figure 3 | Electron energy levels and orbital charge density isosurfaces (isovalue 0.02 coulomb/m 3 ) in 2D CrPc and 2D FePc sheets. The arrows denote 
the spin-up and spin-down. The solid lines correspond to the filled levels and the dotted lines to the unfilled levels. The doubly degenerate levels are 
marked with upper-case letter D. The lines in the rectangles denote levels interacting with 0 2 27i* levels and those in the ovals denote levels interacting 
with CO 5a and 27i* levels. 



of spin conservation 39 " 41 . Nonetheless, for systems such as TMPc 
sheets, the central transition metal atoms exhibit a modest extent 
of spin-orbit effects, enabling them to undergo spin-forbidden tran- 
sitions to achieve a so-called spin inversion. Although the spin -orbit 
interaction is generally regarded to play a key role in mediating the 
spin inversion in reactions containing TM elements, for the 3d trans- 
ition metals, this interaction alone is too small 42 . Therefore, the 
unpaired electrons in 3d channels are predominantly responsible 
for interacting with spin-polarized molecules. For 0 2 molecule, 
which is triplet in its ground state (S = 3), when interacting with 
the CrPc (S = 5) and FePc (S = 3) sheets, the adsorption process 
leads the corresponding adducts to a lower spin state for both 2D 
CrPc-0 2 (S = 3) and 2D FePc-0 2 (S = 1), respectively. Only by 
donating spin up electrons in the 3d channels of Cr and Fe atoms 
to the unoccupied spin down LUMOs (2n*) in 0 2 can such a lower 
spin state be achieved. For a CrPc sheet, the HOMO to HOMO-2 
levels at the T point are all in the a (spin up) state. The doubly 
degenerate HOMO levels predominantly have a d^ character, 
whereas the HOMO-1 and HOMO-2 levels are mainly consisted of 
d z2 and d w components, respectively. All of the four levels lie in a 
notably higher energy range than their counterparts in a 2D FePc 
sheet. The spin up levels in 2D FePc are HOMO-2, HOMO-3, 
HOMO-4 and HOMO-6, as illustrated in Figure 3. In particular, 
the HOMO-3 orbital of a 2D CrPc sheet has a maximized overlap 
with the LUMO (271;*) of 0 2 when the 0 2 molecule adopts a parallel 
(side-on) adsorption configuration due to the symmetry reasons; 
while in a 2D FePc sheet, the dominating interacting orbitals with 
the unoccupied 0 2 2tt* orbitals are HOMO-2 and HOMO-3 exhib- 
iting a d n character, owing to their small energy difference with the 
LUMO of 0 2 . The adsorbed 0 2 molecule will thus adopt an end-on 
configuration with a bent structure where the Fe-O-O angle is max- 
imized for better orbital mixing. As a result of sufficient orbital 
overlapping, more charges will be transferred from the CrPc sheet 
to the adsorbed 0 2 with respect to the case in 2D FePc (Table 1), 
giving rise to an elongated O-O bond length of 1.396 A in a 2D CrPc 
sheet as compared with the value of 1.276 A in a 2D FePc sheet. 



However, the situation completely changes for CO adsorption. 
The central TM atoms interact with the adsorbed CO molecules 
via a well-established electron accepting and electron back-donating 
mechanism to form <j-n coordinate covalent bonds. The <j-n coord- 
inate bond concerns essentially with the donation of CO lone pair 
(5a orbital) to the empty d z2 orbital as well as the back-donation 
from the electron sufficient d n orbital to the electron deficient 2n* 
orbital of CO. Therefore, the energy differences between these per- 
tinent orbitals determine the interaction strength between the CO 
molecule and the TM central atom. In the case of a CrPc sheet, the 
electron-rich d^ orbitals are the doubly degenerate HOMOs lying 
higher in energy, which are facile to donate their electrons to the 
empty CO 271* antibonding orbital. However, the electron-deficient 
d z2 orbital is LUMO + 9, and its energy is too high to achieve 
adequate hybridization with CO 5 a orbital. As a result, the unilateral 
electron donation from Cr to CO results in a negatively- charged 
adsorbed CO, as illustrated in Table 1. On the other hand, the elec- 
tron-rich d K orbitals of a FePc sheet not only contain the spin-up 
orbitals of HOMO-2 and HOMO-3, but also the HOMOs with spin- 
down. Additionally, contrary to the high-lying d z2 orbital of a CrPc 
sheet, the electron-deficient d z2 orbital of the FePc sheet is LUMO + 
1 with a substantially lower energy, thereby enabling CO to inject 
more electrons to it. The bilateral electron injection between CO and 
Fe is responsible for the positively charging of the CO molecule 
adsorbed on a FePc sheet, as listed in Table 1. In summary, the 
pronounced distinctions between the CrPc and FePc sheets in the 
hybridizations between particular d orbitals and the frontier orbitals 
of CO and 0 2 give rise to the opposite adsorption trends for CO and 



Table 1 Hirshfeld charges of CO 


and O2 on 2D CrPc, MnPc and 


FePc sheets 








CrPc 


MnPc FePc 


CO 


-0.0396 


-0.0007 0.0398 


o 2 


-0.3094 


-0.1662 -0.0913 
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0 2 . In the case of CO and 0 2 adsorption on a MnPc sheet, the orbital 
arrangement of a 2D MnPc sheet is intermediate between the CrPc 
and FePc sheets. Therefore, the adsorption energies of CO and 0 2 
also lie between those on 2D CrPc and 2D FePc sheets. Based on the 
above analysis, we select a CrPc sheet for CO oxidation as discussed 
in the following. 

Catalytic mechanism. Surface catalytic reactions proceed via a 
number of elementary adsorption/desorption steps whose reaction 
energetics are associated with the catalytic activity. On the basis of the 
Bronsted-Evans-Polanyi (BEP) relation, the activation energy barrier 
of an elementary reaction step AE can be interpreted as a linear 
function of the dissociative adsorption energy E ad of each adsorb - 
ate involved in the step, providing us a preliminary empirical 
evaluation of the reaction facility 43 ' 44 . Because of the notable E ad of 
0 2 and CO on 2D Pc sheets in conjunction with the superb resistance 
to the CO poisoning effect, we only presented the energy barrier AE 
of the elementary reactions constituting the whole CO oxidation 
cycle on a 2D CrPc sheet in the following part, in order to deter- 
mine the corresponding minimum energy pathway (MEP). We note 
that generally there are two well-established reaction mechanisms for 
the initiation steps involved in surface-catalyzed CO oxidation, 
namely the Langmuir-Hinshelwood (LH) mechanism and the 
Eley-Rideal (ER) mechanism 45 . In the former mechanism, the 
surface co- adsorbed and activated CO and 0 2 molecules interact 
with each other via a concerted reaction pathway to form a peroxo- 
type like intermediate; while in the latter one, the free-standing CO 
directly reacts with a surface activated 0 2 or atomic O generated 
through a surface motivated 0 2 dissociation. Because the 0 2 is 
strongly adsorbed and the O-O bond length is considerably 
elongated upon its adsorption on a 2D CrPc sheet, it seems that 
the ER mechanism would be more favorable. However, contrary to 
our expectation, after full TS search and geometry optimization, we 
find that no C0 2 is formed or other intermediate formation process 
exists, indicating that for this system the CO oxidation is not via the 
ER mechanism. Therefore, the LH mechanism CO + 0 2 — > OOCO 
C0 2 + O followed by the ER mechanism CO + O C0 2 is 
considered. 

Figure 4 and Figure 5 show the energy profiles of MEPs via LH and 
ER mechanisms as well as the side and top views of the local config- 
urations of the adsorbates along the MEPs, respectively. Upon coad- 




-4.52 eV 
FS 



Figure 4 | Schematic energy profile and the corresponding local 
configurations along the Langmuir-Hinshelwood step. Both top (upper 
panel) and side views (lower panel) are illustrated. All energies are given 
with respect to the reference energy of the coadsorbed state for CO and 0 2 
on a 2D CrPc sheet. 



sorption of CO and 0 2 , the CO and 0 2 molecules approach each 
other to form the first transition state (TS1), with an energy barrier of 
0.55 eV along the reaction pathway. Passing over TS1, an peroxide- 
type intermediate state (MS) is formed, which is a slight endothermic 
process with the reaction energy of 0.05 eV for CO + 0 2 — > OOCO. 
Mediated by the central Cr atom, the O-O bond length of the per- 
oxide-type intermediate elongates continuously and passes over the 
second transition state (TS2) with a substantially low energy barrier 
of only 0.09 eV to come to the final state (FS), releasing a large energy 
of 4.52 eV in the overall LH process. The resulting C0 2 molecule is 
easy to desorb from the substrate due to the weaker interaction 
between C0 2 and the substrate, leaving an atomic O strongly inter- 
acting with the Cr atom. Such an activated surface O atom can readily 
bind another CO molecule to form another intermediate state, sur- 
mounting an energy barrier of 0.46 eV. The intermediate state will 
soon turn to the final state without any energy barrier, regenerating a 
2D CrPc sheet with another C0 2 released. The highest energy barrier 
of this catalytic cycle is therefore only 0.55 eV, which is considerably 
lower than that of Fe-N4 porphyrin-like graphene and carbon nano- 
tube at the same calculation level of theory 16 . Such disparate energy 
barriers between Cr and Fe containing porphyrin-like frameworks 
further corroborate the fact that apposite selection of central metal 
atoms can achieve an outstanding performance in catalytic processes. 

Temperature effect. Next we take the temperature effect into 
account to assess the reaction thermodynamics of the induced CO 
oxidation on a 2D CrPc sheet. It is well known that free energy 
change is determined by AG = AH — TAS, where AG denotes the 
Gibbs free energy change in each step throughout the LH and ER 
mechanism, AH is the change in enthalpy, and AS is the change in 
entropy. AH contains two parts, namely A U and APV. In our 
calculation, A U is defined as AU = A£ tot + A£ vi b + A£ rot + AE tram , 
and AS is defined as AS = AS vi b + AS rot + AS tra ns- A£ tot is the 
standard DFT calculated electronic energy at 0 K, and the 
subscript vib, rot and trans are for the contribution of vibration, 
rotation and translation mode, respectively. 

To get further insight into the temperature effect on AH and AS, 
we employ the following equations 




-1.15 eV 
MS 



-1.78 eV 
FS 

Figure 5 | Schematic energy profile and the corresponding local 
configurations along the Eley-Rideal step. Both top (upper panel) and side 
views (lower panel) are illustrated. All energies are given with respect to the 
reference energy of the coadsorbed state for CO and O on a 2D CrPc sheet. 
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Table 2 Free energy changes (in eV) for each step 


in the catalytic CO 


oxidation reaction cycle un 


der different temperatu 


res 


Temperature 




LH Step 




ER Step 




IS-TSl(E al ) 


MS-TS2(E a2 ) 


IS-FS(E rl ) 


IS-TS(E a3 ) 


IS-FS(E r2 ) 


A/ D IS. 


u.oo 


U. 1 z 


A RA 
— 4.00 


n aa 

U.40 


1 PO 

— i .oy 


ZYO. 1 J IX 


n ra 
u.oo 


U. 1 z 


A RA 
— 4.00 


n aa 
U.40 


1 o 

— i .y 


300 K 


0.56 


U. 1 z 


A RA 
— 4.00 


n /1A 
U.40 


1 o 

— i .y 


325 K 


0.57 


U. 1 o 


A RA 
— 4.00 


n /ia 
U.40 


1 O 1 

— i .y i 


350 K 


0.57 


0.13 


-4.56 


0.46 


-1.92 


375 K 


0.57 


0.13 


-4.56 


0.46 


-1.93 


400 K 


0.58 


0.14 


-4.57 


0.46 


-1.94 


425 K 


0.58 


0.14 


-4.57 


0.47 


-1.95 


450 K 


0.58 


0.15 


-4.57 


0.47 


-1.96 


475 K 


0.59 


0.15 


-4.57 


0.47 


-1.97 



As(r 2 )=As(r,)+ 



,Ev B c p (Byr 



T 2 



(2) 



where AC p denotes the change of heat capacity at constant pressure 
and v B denotes stoichiometric number. The corrected free energy 
changes for the whole CO oxidation on a 2D CrPc sheet under 
different temperatures (275 K ~ 475 K) are listed in Table 2. We 
can see that although elevating the temperature can make the whole 
reaction thermodynamically more viable, the increased energy 
barriers for surmounting the three transition states hinders the 



dynamical feasibility. The activation energy barriers alter only 
slightly in the temperature range between 275 ~ 475 K, suggesting 
the notable stability and catalytic feasibility at lower temperatures as 
well as higher temperatures. 

Electronic structure analysis. To unveil the microscopic scenario 
behind the peculiar catalytic activity of the 2D CrPc sheet, we 
calculated the local density of states projected onto the 01-02 and 
the C-O bonds as well as the d-projected LDOS of Cr in IS, TS1, MS 
and TS2 of the LH step, respectively. As illustrated in Figure 6, the 
271* level of 0 2 is partially shifted downward upon the adsorption in 
the IS state due to the notable charge transfer from the Cr atom. In 




-12 -10 -8 -6 -4 -2 0 2 4 6 -12 -10 -8 -6 -4 -2 0 2 4 6 



Energy (eV) Energy (eV) 




-12 -10 -8 -6 -4 -2 0 2 4 6 -12 -10 -8 -6 -4 -2 0 2 4 6 



Energy (eV) Energy (eV) 

Figure 6 | Spin-polarized local density of states (LDOS) projected onto C-O and Ol -02 on a 2D CrPc sheet, as well as the d-projected LDOS of Cr in the 
IS, TS1, MS and TS2 in the LH step. Black solid curve for the LDOS projected on C-O; blue solid curve for the LDOS projected on 01-02; red solid curve 
for the d-projected LDOS of the Cr atom. The corresponding orbital notations are also assigned. 
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contrast, the 2n* level of CO is still empty because of the paltry 
amount of charge transfer between CO and Cr. For the 01-02 
species on the 2D CrPc sheet, from IS to TS1, owing to the 
electronic resonance with CO 27r* state, the 2 71* level of 0 2 is 
further populated and broadened to hybridize with Cr 3d state to a 
larger extent. Meanwhile, the 2n* level of CO is also partially 
populated in the TS1 state. The 01-02 In state is elevated and 
becomes more involved with the Cr 3d state from TS1 to MS to 
TS2, which is ancillary in weakening the 01-02 bond. The 271* 
level of 01-02 is gradually elevated upward as a result of the 
generation of the new C-02 bond from TS1 to TS2. For the C-O 
species on the 2D CrPc, because of the end-on adsorption 
configuration, the 5 a level is shifted upward and strongly 
hybridized with the Cr 3d level in the whole LH step rather than 
the In level, which is one of the predominant factor in the 
dissociation of the 01-02 bond because of the side-on adsorption 
configuration of 0 2 molecule. Furthermore, the magnetism of the 
entire system is quenched at the TS2 state originating from the 
formation of the nonmagnetic C0 2 molecule and O adsorbed 2D 
CrPc sheet. Overall, mediated by the Cr 3d state, the 01-02 In and 
271* level interact strongly with the C-O 5a and 2tt* level. These 
electronic state resonances lead to the formation of C-02 bond 
and the dissociation of 01-02 bond, accelerating the catalytic 
elimination of the hazardous CO molecules. 

Discussion 

Based on first-principles calculations, we find that the 2D CrPc sheet 
exhibits superior catalytic activity towards CO oxidation withstand- 
ing CO poisoning effect. The exceptional catalytic activity can be 
ascribed to its peculiar orbital order, which enables the 2D CrPc sheet 
to interact strongly with 0 2 while adsorbing CO only to a modest 
extent. The spin-up occupied d^, d xy and d z2 orbitals are majorly 
accountable to the electron injection to the 0 2 void 2n* orbital. 
On the other hand, only three occupied d^ orbitals provide electrons 
to partially fill the empty CO 27t* orbital, the unoccupied d z2 orbital 
lies in a too high energy range to interact effectively with the CO 5 a 
orbital. Furthermore, we identified that the CO oxidation reac- 
tion catalyzed by the CrPc sheet proceeds first via Langmuir- 
Hinshelwood (LH) mechanism and then via Eley-Rideal (ER) 
mechanism with highest energy barrier of 0.55 eV, much lower than 
its Fe- containing counterparts, suggesting this reaction is viable to 
proceed at low temperature. Moreover, our vibrational analyses reveal 
that the 2D CrPc framework is viable in catalyzing the whole CO 
oxidation reaction at a wide temperature range without being signifi- 
cantly affected by temperature variation. The high activity of the 2D 
CrPc sheet along the LH step may be attributed to the electronic 
resonance among electronic states of CO, 0 2 , and the Cr atoms, 
particularly among the Cr-3d, CO-5a, 00-271*, and 02-2tt* orbitals. 

For avoiding the CO poisoning effect, PtAu alloy bimetallic cata- 
lysts and Pt overlayer-metal carbide catalysts have been suggested to 
improve the catalytic efficiency and lower the adsorption energy of 
CO to a moderate extent 46 ' 47 . However, the low surface to volume 
ratio impedes further improvement to their catalytic activity. 
Although atoms or clusters loaded on pristine or defective carbon 
based material such as graphene showed high surface to volume 
ratios, previous calculations suggested that almost all of them suffer 
from different extent of CO poisoning effect. Due to the following 
advantages, the studied 2D CrPc sheet is superior to platinum based 
catalysts and carbon based nanocomposite catalysts: (1). The syn- 
thetic procedure of 2D TMPc sheets includes polymerization of small 
TCNB molecules and TM atoms, which creates no byproducts and is 
easily transferable to both metallic and insulating substrates, which is 
highly desirable in realistic device fabrication; (2). It possesses cata- 
lytically active unsaturated metal ions with lower valence states, and 
it exhibits comparable amount of available active metal sites with that 
in metal nanoparticles. To compare the catalytic efficiency of the 2D 



CrPc sheet with nanoparticles, we use the available catalytically 
active metal site number per gram instead of surface to volume ratio. 
In 2D CrPc, each structural unit contains 20 carbon atoms, 8 nitro- 
gen atoms, 4 hydrogen atoms and 1 chromium atom. Therefore, The 
number of active metal sites in 1 gram of 2D CrPc sheet is 

Ni = N A = x 6.02 x 10 23 = 1.48 x 10 21 . Which is compar- 

408 408 Y 
able to the numbers of 1.87 X 10 21 and 1.78 X 10 21 for cuboctahedral 
Pt 20 i and icosahedral Au 55 nanoparticles, respectively. (3). Cr atoms 
are uniformly distributed in 2D CrPc sheet without the problem of 
clustering which is often encountered in other substrates loaded with 
catalysts; (4). Different from the broad d band features of platinum 
based catalysts with nanocrystal structure, the d band of the sparsely 
distributed Cr atoms is sharp near the Fermi level and can be effec- 
tively tuned by other ways such as strains or chemical modification. 
These features discussed above make the 2D porous CrPc sheet very 
unique and promising in catalytic CO oxidation for large scale appli- 
cations in industry. 

Methods 

Periodic spin-unrestricted density functional theory calculations were performed 
using the DMoP code 48 ' 49 . The generalized gradient approximation (GGA) with the 
Perdew — Burke — Ernzerhof (PBE) functional was employed to describe the 
exchange and correlation potential 50 . DFT semicore pseudopotentials (DSPPs) 51 were 
used, where all-electron calculations were performed for C, O and N atoms, while the 
core electrons were replaced by a single effective pseudopotential including the 
relativistic effects for TM atoms. Double numerical plus polarization (DNP) basis 
set 48 was employed and the real space global orbital cutoff radius was set as 4.8 A to 
ensure high quality results. We used a PBE-D scheme with the Grimme vdW cor- 
rection 5253 , which has been well validated to describe the adsorption of small mole- 
cules on 2D nanomaterials 54 . 

We applied periodic boundary conditions with a vacuum space of 15 A to avoid 
interactions between the TMPc sheet and its adjacent periodic images. A single unit 
cell was used to model the 2D TMPc systems because of the negligible effects of 
magnetic coupling between the unit cells on their catalytic properties. For geometry 
optimization and the search for the transition state (TS), the Brillouin zone was 
sampled with 5 X 5 X 1 k points. For the electronic structure calculations, a 
Monkhorst- Pack scheme with 19 X 19 X 1 k points grid meshes was employed. The 
atoms were fully relaxed without any symmetry constraints, with the convergence in 
energy, force and displacement being set at 1 X 10" 5 Ha, 0.002 Ha/A and 0.005 A, 
respectively. All the self- consistent field (SCF) calculations were performed with the 
convergence criterion of 1 X 10~ 6 Ha without any smearing. The transition states 
were located through a complete LST/QST protocol involving linear synchronous 
transit (LST) maximization followed by repeated conjugated gradient (CG) refine- 
ments, and then several quadratic synchronous transit (QST) maximizations with 
repeated conjugated gradient (CG) refinements between them 55 . TS optimizations 
were carried out to confirm that the located TSs contain only one imaginary fre- 
quency. The nudged elastic band (NEB) method 56 ' 57 was used to verify the MEP for 
elementary reaction steps. 
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